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Abstract 

I argue that particles heated by relativistic shocks should assume an equilibrium energy 
distribution. This leads to a synchrotron spectrum F u <x i/ 1 / 3 up to approximately the 
critical frequency u 0 of an electron with the mean electron energy. Application to GRBs 
implies that a burst with 10~ 5 erg/cm 2 s of soft gamma-rays and hu 0 = 300 KeV should 
be about 18th magnitude in visible light and a few /ijy at 1 GHz (less if self-absorbed). 
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1. Introduction 


t 


A model of gamma-ray bursts (GRB) has been proposed in which a relativistic fire- 
ball’s debris shell (Shemi and Piran 1990) forms a relativistic collisionless shock when 
it interacts with surrounding low-density gas (Rees and Meszaros 1992, Meszaros and 
Rees 1993a, Katz 1994a). The shock-heated plasma radiates by the synchrotron process, 
producing the observed gamma-rays. It has been usual in astrophysics to assume that 
shock-accelerated particles have a power-law distribution in energy. In this paper I argue 
that this assumption is inapplicable to relativistic shocks, for which I predict a thermal 
equilibrium particle distribution function and a synchrotron spectrum F„ oc i/ 1 / 3 below its 
peak. 


2. Fermi Acceleration 

Most particle acceleration mechanisms considered in astrophysics (with the exception 
of acceleration in a large quasi-static electric field in pulsars) axe based on developments 
of Fermi’s (1949) model of stochastic acceleration. These mechanisms include acceleration 
by a nonrelativistic shock (Blandford and Eichler 1987), turbulent plasma acceleration 
(Sagdeev and Galeev 1969) and viscous shear acceleration (Katz 1991). In these pro- 
cesses the accelerated particles are not in thermal equilibrium, and the evolution of their 
distribution function is described by a Fokker-Planck equation. 

Fermi derived the general result that the distribution function produced by these 
stochastic acceleration processes is proportional to a power of the particle energy, but a 
dimensional argument is sufficient: Any deviation from a power law would define a char- 
acteristic energy, yet in these models of stochastic acceleration there are no characteristic 
energies between a particle’s rest mass and an extremely high energy above which it is not 
magnetically confined to the region of acceleration. The number density n of accelerated 
particles is indeterminate because there is an essentially infinite reservoir of (nonrelativis- 
tic) thermal particles. The energy density e of accelerated particles is also indeterminate 
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because the acceleration is driven by a reservoir of fluid kinetic energy in the shock or 
plasma or hydrodynamic turbulence, an indeterminate fraction of which is converted to 
particle acceleration. The prediction of a power-law differential distribution of particle 
energy under these conditions 

N(E) oc E~ p , (1) 

where N(E) is the number of particles per unit energy per unit volume, is borne out by 
observations of cosmic rays and of power-law synchrotron spectra radiated by accelerated 
particles 

F u oc u~ 3 (2) 

in many astronomical objects; in optically thin synchrotron theory (Rybicki and Lightman 
1979) s = (p- l)/2 if p > 1 /3. 

In order that the total number and energy of accelerated particles be finite, a distri- 
bution of the form (1) must have at least one characteristic energy E 0 at which p changes, 
with p = p< < 1 for E < Eo and p — p > > 2 for E > Eq . The corresponding synchrotron 
spectral indices s< < 0 and s> > 1/2. In most astronomical synchrotron sources (but not 
GRB) E 0 is very small (the rest mass energy, or the injection energy below which Coulomb 
drag is significant) and outside the range of observation, and the spectral index s< < 0 
is not observable (synchrotron radiation by such low energy particles is negligible or at 
unobservably low frequencies, and the particles may be nonrelativistic, for which different 
results apply). Observed spectral indices usually exceed 1/2 (for cosmic rays the observed 
p > 2), as expected for E > Eq. 


3. Relativistic Shocks 


It is often assumed that a power-law particle distribution function (Eq. 1) will apply 
to any collisionless gas of shock-accelerated particles. This is incorrect for relativistic 
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collisionless shocks*. In such a shock n and e are determined by the jump conditions, and 
the mean energy per particle sets the natural energy scale Eq = e/n. 

In order for a shock to occur there must be an irreversible process, which mixes the 
distribution function in phase space, increasing a coarse-grained entropy. Without detailed 
knowledge of the mechanics of this process, it is natural to assume that it leads to a final 
state in which all microstates consistent with the constraints (on n and e) axe equally prob- 
able. This assumption defines a microcanonical ensemble (Reif 1965), and was also made 
by Lynden-Bell (1967) in his discussion of violent gravitational relaxation. The resulting 
distribution of particle energies, assuming relativistic kinematics and nondegeneracy, is 
that of thermal equilibrium 


N(E) oc E 2 exp(-3E/Eo); (3) 

the “temperature” ksT = E 0 /3, and the parameters p< = —2 and p > — > oo. This 
argument for (3) depends essentially on the existence of the constraints. 

In a colhsionless shock (3) must result from the interaction of coarse-grained clumps 
in phase space, rather than from single-particle collisions. These clumps interact electro- 
magnetically through plasma turbulence, and we may regard the fields as the means by 
which the clumps interact, in analogy with the gravitational interaction of coarse-grained 
clumps considered by Lynden-Bell. 

If the phase space clumps are sufficiently long-lived (as, for example, may be the 
clumps discussed by Dupree [1982]), they may come to equilibrium with each other and 
with the turbulent fields, considered as independent degrees of freedom. This hypothesis 
is distinct from that of violent relaxation. The brightness temperature of the turbulence 

* The argument presented by Katz (1994a) for p = 3/2 and s = 1/4 is wrong; nonrel- 
ativistic shock acceleration theory is inapplicable, and the assumed compression ratio is 
inconsistent with the (correct) values cited elsewhere in that paper. 
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is given by 


k B T> = (4) 

where T v is the spectral density of the turbulence. In a relativistic plasma most modes 
will have relativistic phase and group velocities, so we may approximate 




B 2 c 

87rAi/’ 


( 5 ) 


where B is the turbulent magnetic field and Ay is its spectral bandwidth. A clump of size 
A will contain « n\ 3 particles and a total energy ~ u\ 3 Eq « ksT c , where T c is the kinetic 
temperature of the clumps. Taking Ai/«i/ and A « c/v and equating T c and Tj yields 

B 2 

e = nE 0 « — . (6) 


This demonstrates that equipartition between clumps and plasma turbulence is consistent 
with the energetics, and lends credibility to the equipartition of particle and magnetic 
energy assumed by Katz (1994a). 

The length A is determined by the wavelengths of the fastest growing plasma insta- 
bilities. These may be two-stream instabilities when the plasmas first interpenetrate, and 
then velocity space anisotropy instabilities. Equation (6) also implies approximate equal- 
ity between the plasma and gyro-frequencies u> p and u> g , so that the distribution of A is 
probably peaked around A « 2irc/u> p ~ 27rc/u> 5 , with T v peaked around v « c/A. 

Electron-ion equipartition cannot result from the very slow single-particle interactions, 
and must result from electromagnetic coupling between charge- or current-unneutralized 
clumps in electron and ion phase space. These complex processes are beyond the scope of 
this paper. 

Finally, a nonequilibrium distribution can also be shown to relax to the equilibrium (3) 
if an H-theorem is satisfied, which follows from the assumption of molecular chaos (Liboff 
1969). This assumption holds if the phase space clumps interact weakly, but is likely to be 
more general — H-theorems are observed to hold in dense and strongly coupled systems for 
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which it is hard to justify the assumption of molecular chaos. Strong two-particle (or two- 
clump) momentum correlations axe disrupted, at least in ensemble average, as the particles 
(or clumps) interact with many others, and thermodynamic equilibrium is observed even 
though the assumption of molecular chaos is hard to defend. 

4. Predicted Spectra 

For electron distributions with p < 1/3 the synchrotron spectrum at frequencies below 
the critical frequency is dominated by the radiation of the most energetic electrons (Jackson 
1975), and s = —1/3 (rather than (p — l)/2). Thus the particle distribution (3) leads to a 
predicted radiation spectrum 

F v oc i/ 1 / 3 . (7) 

This power law spectrum survives averaging over emission regions with a range of Eo, 
magnetic field, and Doppler shift if the frequency of observation is everywhere below the 
critical synchrotron frequency (Doppler-shifted to the observer’s frame) for electrons with 
the local Eo in the local field. The observation (Schaefer 1994) in GRB at soft X-ray 
energies of N„ oc u~ 07 , corresponding to F„ oc u 03 , may be confirmation of (7). The 
steeper spectra of GRB at harder X-ray and gamina-ray energies may reflect evolving 
heterogeneous source regions with a distribution of (Doppler-shifted) critical synchrotron 
frequencies, some of which are below the frequency of observation. 

In GRB the spectrum (7) should extend from X-rays to radio frequencies, where self- 
absorption becomes significant (Meszaros and Rees 1993b, Paczynski and Rhoads 1993, 
Katz 1994ab). Because the radiation at all frequencies is dominated by electrons with 
E « Eo, in the self-absorbed region s = — 2 rather than the value s = —5/2 found 
(Rybicki and Lightman 1979) when p > 1/3. 

The spectrum of GRB may be normalized to the observed soft gamma-ray fluxes. The 
normalization is necessarily very approximate, because source regions are likely heteroge- 
neous and the transition between (7) and the expected exponential cutoff at high photon 
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energies may therefore extend over several decades of the integrated spectrum, including 
the soft gamma-ray region. An intense GRB with a flux 10 -5 erg/cm 2 s in a bandwidth of 
400 KeV around hu = 300 KeV has a flux there of « 10 mJy. Extrapolation leads to a 
flux of « 0.2 mJy in visible light, (an ~ 18th magnitude star) and to w 2 //Jy at 1 GHz. 
The effective flux of a brief transient measured by a broad-band receiver may be further 
reduced by plasma dispersion. Because these fluxes are low (and, in addition, self-absorbed 
at low frequencies) intergalactic dispersion (Ginzburg 1973, Palmer 1993, Katz 1994ab) 
and optical counterparts to GRB will be difficult to observe. For an actual GRB it is more 
accurate to extrapolate downward from the highest frequency vq at which (7) is observed, 
which may lead to substantially higher visible and radio fluxes if s > —1/3 between vq and 
soft gamma-ray frequencies. 

After the initial gamma- ray transient, a relativistic fireball will continue to expand and 
to produce radiation up to a critical frequency which decreases with time as the blast wave 
degrades (Paczynski and Rhoads 1993, Katz 1994ab, Meszaros and Rees 1993b, 1994). 
In a uniform medium (a very unrealistic assumption, as shown by the complex temporal 
structure of observed GRB) the intensity at a given frequency will increase with time (in 
a simple model oc < 4 ' 5 ; Katz 1994b) until the critical synchrotron frequency decreases to 
the frequency of observation. However, at all times the spectrum (7) should be observed 
in the very broad range of frequencies between the critical frequency and the onset of 
self-absorption. 


5. Discussion 

This paper makes two predictions for the characteristic spectrum of synchrotron radi- 
ation produced by relativistic shocks. The first prediction, s< < 0, uses only the finiteness 
of the total particle number. The second, more specific, prediction, s < = —1/3 (7), is 
based upon the argument (§3) for an equilibrium particle distribution function (3), though 
any distribution with p < 1/3 is sufficient to lead to s — —1/3. 
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The equilibrium particle distribution (3) also predicts an exponential cutoff on the ra- 
diation spectrum at high frequency. This has not generally been observed in GRB, which 
may be explained if the observed radiation is the superposition of emission from a distri- 
bution of emitting regions with a wider range of parameters and characteristic synchrotron 
energies, extending up to the the highest photon energies observed. In addition, because 
of the poor counting statistics at high photon energies, particularly when there is a rapid 
cutoff, it is hard to verify the functional form of the cutoff. 

The arguments presented here for relativistic shocks are general, and not specific to 
GRB. For example, they should apply to relativistic blast wave models of AGN such as 
those proposed by Blandford and McKee (1977). 

I thank P. Diamond and B. E. Schaefer for discussions and NASA NAGW-2918 for 
support. 
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